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i^ I The recent theoretical discovery of finite-aniphtude travelhng waves in pipe 

^^ • flow has re-ignited interest in the transitional phenomena that Osborne Reynolds 

^ I studied 125 years ago. Despite all being unstable, these waves are providing fresh 

^^ ' insight into the flow dynamics. Here we describe two new classes of /ii(7/i??/-symnietric 

^ ' I travelling waves (possessing rotational, shift-&-reflect and mirror symmetries) and 

report a new family of mirror-symmetric waves which is the first found in pipe 

M i flow not to have shift-&-reflect symmetry. The highly-symmetric waves appear at 

>^p>' lower Reynolds numbers than the originally-discovered non-mirror-symmetric waves 

found by Faisst & Eckhardt 2003 and Wedin & Kerswell 2004 and have much higher 

^ ■ wall shear stresses. The first M-class comprises of the various discrete-rotationally- 

t^_| I symmetric analogues of the mirror-symmetric wave found in Pringle & Kerswell 

5/3 ■ (2007) and have a distinctive double layer structure of fast and slow streaks across 

O I the pipe radius. The second N- class has the more familiar separation of fast streaks 

C/3 ■ to the exterior and slow streaks to the interior and looks the precursor to the class 

^^1 of non-mirror-symmetric waves already known. 

f^ ' Keywords: Pipe flow; travelling waves; turbulence transition 

1. Introduction 
> , 
'sj" I The stability of the flow of a Newtonian fluid such as water through a straight pipe of 

"^ ' constant circular cross-section has fascinated scientists ever since Reynolds' (1883) 

original experiments. Reynolds realised that there was a single non-dimensional 
control parameter for the flow, Re := UD/v (where U is the mean velocity, D is 
the pipe diameter, and v is the kinematic viscosity of the fluid), but that there 
was no unique value of this parameter for transition to occur. Instead, he noticed 

^**^ , that the value of Re required for transition depended critically on how cleanly he 

performed the experiment, that is, on the ambient level of noise. Later theoreti- 
cal and computational work confirmed that the steady, unidirectional, parabolic, 
laminar 'Hagen-Poiseuille' flow (Hagen 1839, Poiseuille 1840) realised at low Re is 
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H \ linearly stable at least to Re — 10^ (Meseguer & Trefethen 2003) and that pipe 

flow transition must be a finite-amplitude process. The variability in the critical Re 
across experimental realisations has indicated that this process is also sensitive to 
the exact form of disturbances present in the flow (e.g. Peixinho & MuUin 2007). 

Once triggered, transition is generally abrupt and can lead to both localised and 
global forms of turbulence depending on Re. For 1760 < i?e < 2, 300, turbulence 
is localised into 'puffs' (Wygnanski & Champagne 1973) which have a weak front, 
a sharp trailing edge, and maintain their length at about 20 D as they propagate 
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along the pipe. For 2, 300 < Re, the turbulence delocalises into 'slug' turbulence 
(Wygnanski et al 1975) which has sharp front travelling faster than the mean flow 
and a trailing edge travelling slower so that it expands while propagating down- 
stream. As a result of this varied behaviour and the sensitivity of the transition 
process, pipe flow has become the canonical example of shear flow transition. 

From a dynamical systems viewpoint, the existence of turbulence requires the 
presence of simpler solutions of the governing equations to provide a framework in 
state space to sustain turbulent trajectories from simply relaminarising. Natural 
questions are then whether the emergence of such alternative solutions can be used 
to predict the threshold (lowest) Re for transition and what role they play in the 
process as well as the final turbulent state. Steady states and travelling waves (TWs) 
- fixed points in an appropriately chosen Galilean frame - are the simplest solutions. 

The first step towards identifying such states was taken by Smith & Bodonyi 
(1982) who predicted the existence of helical modes (with swirl) through a critical 
layer analysis. Landman (1990), however, failed to find any numerical evidence for 
their existence. In the same year, Nagata (1990) used a homotopy approach to find 
finite amplitude solutions in plane Couette flow starting from the linearly-unstable 
rotating situation. A similar attempt using rotating pipe flow, however, failed to 
reach the vanishing rotation limit (Barnes & Kerswell 2000). By this time a more 
systematic approach was emerging based upon the idea of a self-sustaining process. 
Hamilton et al. (1995) suggested that streamwise rolls could create streamwise 
streaks which in turn could reinforce the original rolls by some unspecified non- 
linear process. Waleffe (1997,1998) closed the loop by showing that the streaks had 
infiectional instabilities which led to what he described as 'wriggling' in such a man- 
ner that energy is fed back into the initial rolls. He then identified such solutions 
in plane-Couette flow using a constructive homotopic approach. Faisst & Eckhardt 
(2003) and Wedin & Kerswell (2004) employed this technique within the geometry 
of pipe flow to flnd the expected travelling waves although all were without swirl 
and therefore unrelated to the predictions of Smith & Bodonyi (1982). Experimen- 
tal observations (Hof ct al. 2004,2005) and numerical computations (Kerswell & 
Tutty 2007, Schneider et al. 2007a, Willis & KersweU 2008) then followed which 
indicated that these waves are realised albeit transiently as coherent structures 
in turbulent pipe flow. These waves are all born in saddle node bifurcations with 
the lowest at Re = 1251 for TWs with a discrete rotational symmetry. Pringle & 
Kerswell (2007) substantially lowered this to 773 by finding travelling waves with 
no rotational symmetry. The Re gap between when alternative solutions exist and 
when sustained turbulence occurs is then certainly positive but also too large to be 
a useful predictor. 

Recently Duguet et al (2008a) have adopted an entirely different approach to 
finding TWs based upon identifying episodes in time-dependent pipe fiow which 
exhibit temporal periodicity. To generate a long time signal for this search, sus- 
tained fiow dynamics away from the laminar state is required. The obvious source 
for this is studying flow on the turbulent attractor but they realised that there was 
a more promising alternative where the flow is less energetic and periodic episodes 
are more clearly evident: the laminar-turbulent boundary. This is the set of flows 
which separates initial conditions which immediately relaminarise from those which 
undergo a turbulent episode. An initial condition chosen on this boundary or 'edge' 
(Skufca et al. 2006) will evolve in time forever treading a tight rope between relam- 
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inarising and becoming turbulent (Toh & Itano 1999, Itano & Toh 2001, Schneider 
et al 2007b). Duguet et al (2008a) identified 3 new TWs by using velocity fields 
from near-time periodic episodes in this energetically intermediate flow dynamics 
as starting guesses for an iterative Newton-Krylov solver. 

The purpose of this paper is to show how the TWs found by Duguet et al (2008a) 
led to the discovery of 2 new classes of mirror-symmetric TWs hereafter referred to 
as the M- and N-classes, and a distinctly new family of mirror-symmetric TWs. Each 
class is partitioned by the discrete rotational symmetry (see (2.1) below) enjoyed by 
the waves which defines a family within the class. This suggests a natural labelling 
system where, for example, the family of M-class TWs with m (integer) degree of 
discrete rotational symmetry will be referred to as the Mm family. We argue that 
the mirror-symmetric wave of Pringle & Kerswell (2007) is actually part of the first 
family (Ml) of the M-class and 'A3' in the nomenclature of Duguet et al (2008a) 
is a member of the second family (M2). We further report the families M3, M4, 
M5 and M6 which, unlike anything seen so far, all have a double layer structure of 
fast and slow streaks across the pipe radius. The second N-class of TWs has 'C3' in 
Duguet et al (2008a) as part of its second family (N2) with the families N3, N4 and 
N5 described here for the first time. In these waves, the fast streaks are positioned 
near the pipe wall and slow streaks in the interior, a separation familiar from the 
original non-mirror symmetric TWs found by Faisst & Eckhardt (2003) and Wedin 
& Kerswell (2004). Since these waves only have shift-&-reflect symmetry - see (2.2) 
below - we refer to this original set as the S-class hereafter for convenience. 

The plan of the paper is to start by discussing the various symmetries of the 
TWs which serve to categorize them in section 2. Beyond these, all the TWs are 
parametrised by their axial wavenumber a which is a continuous variable with a 
finite range at a given Re. These waves are exact solutions of the Navier-Stokes 
equations when periodic flow conditions are imposed across the pipe with a spatial 
period equal to an integer multiple of l-njo.. The 3 TWs found recently by Duguet 
et al (2008a) are then introduced in section 3. Numerical codes developed in Wedin 
& Kerswell (2004) are used to explore these waves in parameter space. A new 
classification is then introduced in section 4 stimulated by the discoveries of related 
families. A final discussion follows in section 5. 



2. Travelling wave symmetries 

The original work by Faisst & Eckhardt (2003) and Wedin & Kerswell (2004) used 
the homotopy approach developed by Waleffe (1997) to look for exact steady solu- 
tions in an appropriately chosen Galilean frame. In this, an artificial body forcing 
is added to the Navier-Stokes equations in order to drive streamwise rolls. These 
advect the underlying shear towards and away from the pipe wall, thereby gener- 
ating positive and negative streamwise anomalies in the mean flow called 'streaks'. 
If these streaks are of sufficient amplitude they become inflexionally unstable to 
streamwise-dependent flows which if allowed to grow to sufficient amplitude can re- 
place the artificial forcing as the energy source for the rolls. By choosing the forcing 
functions to have certain discrete rotational symmetries defined by 

R„i : (m, u, w,p)(s, (/), z) — > (u, u, w,p)(s, (/) + 27r/m, z), (2-1) 
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and selecting streamwise-dependent instabilities having the shift-&-reflect symme- 
try 

S: (u, u, i(;,p)(s, (/), z) ^ (u, — w, w,p)(s, — (/), z + Tr/a), (2-2) 

where lix j a is the wavelength, exact TWs were found by continuing back to zero 
forcing in ?7i = 2,3,4,5 and 6 discrete rotational symmetry subspaces. All appear 
through saddle node bifurcations with the lowest at Re — 1251 corresponding to a 
R3 symmetric TW. Each TW family is parameterised by its axial wavenumber a 
which, at a given _Re, has a finite range (e.g see figures 1 & 2 of Kerswell & Tutty 
2007). 

Later, motivated by a chaotic flow state found on the laminar-turbulent bound- 
ary by Schneider et al (2007b), Pringle and Kerswell (2007) generalised the forcing 
function to have no rotational symmetry (see also Mellibovsky & Meseguer 2007). 
As a result, they found an asymmetric TW (formally the missing Ri state possess- 
ing one pair of streamwise rolls) with S symmetry. This TW originates through a 
symmetry-breaking bifurcation from a state which satisfies the additional shift-&- 
rotate symmetry 

Hjn : (u, u, i(;,p)(s, (/), z) ^ (u, u, w,p)(s, (/) + tt/to, z + tt/q:) (2.3) 

with TO = 1 and is therefore mirror-symmetric about = 7r/2 (the plane = 
being set by the shift-&-reflect symmetry). For general ?Ti, 

SrJin^^Zj^ j = l,3,5,...2m-l (2.4) 

where rim'' implies Hm is applied j times and 

Z^ : (u,v,w,p)(s,0, z) -> (m, -u,u;,p)(s,2V' - (/>, z) (2.5) 

represents reflection in the plane (/) — ijj. Helical generalisations, which satisfy a 
modified shift-&-rotatc symmetry 

n'^ : (u,v,w,p){s,4>,z) -^ {u,v,w,p)(s,(j) + (I ^ )7r, z + tt/q:) (2.6) 

a 

(where f3 is the helicity, Pringle & Kerswell 2007) but no S-symmetry, were also 
found but always at higher Re indicating that the fiow prefers the streaks to be 
aligned with the flow rather than twisted around it. 

3. Travelling waves found on the laminar-turbulent 

boundary 

The new TWs discovered by Duguet et al (2008a) were found by studying the 
flow dynamics on the laminar-turbulent boundary within the R2-symmetric sub- 
space. Surprisingly, at Re = 2875 and periodic pipe length 5 D, the rotationally- 
unrestricted situation only ever revealed the already-known asymmetric TW (Pringle 
& Kerswell 2007) and a weakly rotating version of it (A in Duguet et al 2008a). 
The fact that new R2-symmetric TWs were observed as transient coherent struc- 
tures only in the R2-symmetric subspace calculations is due to the reduced number 
of unstable directions they have as non R2-symmetric flows are removed from the 

Article submitted to Royal Society 



Highly- symmetric travelling waves in pipe flow 




Figure 1. The travelling waves ^3 (top left), C3 (top right) and Z2 (bottom left). For 
comparison purposes we also include 5*2 (26_1.25 in Kerswell & Tutty 2007). All are shown 
at their respective values of a* (see Table 1) and at Re = 2400 except for Z2 which is 
shown at Re = 3860. The contours indicate the magnitude of the streamwise velocity 
difference from the underlying laminar flow. The contours and colouring is standardised 
across the plots (and more generally figures 3, 4 and 7) with contours increments of 0.02(7 
running from — 0.19t/(dark red) to 0.19C/(white) (the colour outside of the pipes indicates 
zero). Arrows indicate the cross-stream velocities. 



dynamics. The chance of the flow 'visiting' a TW in phase space is presumably re- 
lated to how unstable it is (the spatial periodicity of the flow was also shortened to 
improve the stability of the TWs). In the notation of Duguet et al (2008a) the 3 new 
TWs were labelled as 'A3', 'C3' and 'D2' and all are mirror symmetric even though 
this symmetry was not imposed on the flow. A3 and C3 were readily converged in 
the continuation codes of Wedin & Kerswell (2004) but D2 could not be, due, we 
suspect, to insufficient axial resolution being achievable in the continuation code. 
The roll structure of D2 was, however, used to design a forcing in the homotopy 
approach which successfully yielded a similar-looking TW - called Z2 - with the 
same symmetries as D2. That they were, in fact, different waves became apparent 
when Z2 could not be continued below Re = 3250 whereas D2 was discovered at 
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Re = 2875. The structure of A3, C3 and Z2 is shown in figure 1 compared to a 
known TW in S2, which, at a = 1.25 and Re ~ 2400 is already known to be on 
the laminar-turbulent boundary (this is 26_1.25 in the nomenclature of Kerswell & 
Tutty 2007). 

While all the new families possess the apparently universal features of exact co- 
herent structures known in wall-bounded shear flows - wavy streaks with staggered 
quasistreamwise vortices (Nagata 1990, Waleffe 1998, 2001, 2003, Faisst & Eckhardt 
2003, Wedin & Kerswell 2004, Pringle & Kerswell 2007), there are new structural 
features. A3 has a strikingly different cross-sectional profile from previously known 
TWs (e.g. S2) in that both fast and slow streaks are concentrated at the pipe wall 
leaving the interior relatively quiescent. Figure 1 shows that C3 and Z2 are also 
noticably different too. However, it is the axial structure of Z2 which really makes it 
stand out. Z2 is not S-symmetric nor rj-symmetric but does have mirror symmetry 
about two perpendicular planes. There is no a priori reason to expect the flow to 
prefer one of these symmetries over any of the others, but to date all TWs in any 
of the canonical shear flows have always been S-symmetric (except, trivially, the 
helical modes of PK07; Waleffe, private communication). Therefore Z2 and D2 are 
the first TWs found to possess only Z-symmetry (Z^ with tp suppressed as there is 
no longer an origin for </>). Both Z2 and D2 are, however, close to being S-symmetric 
in the sense that the simple indicator 

L fn (u — Su)^L-osdsd(/) , Q, , , 

■'%^-'y -^^^ - = 0(10-3). (3.1) 

/q''/q (u + Su)2|2^osdsd0 

In contrast, A3 and C3 are both S- and J72-symmetric and hence also reflec- 
tionally symmetric about a plane at ±7r/4 to the plane of shift-&-refiect symmetry 
(horizontal in figure 1) i.e. they have Z-|-7r/4 symmetry. The modes C3 and A3 appear 
through saddle node bifurcations at Re — 1141 and 1125 with optimal wavenum- 
bers a* = 1.2 and 2.0 respectively giving these lowest saddle node bifurcations: see 
Table 1. 

(a) Stability 

An interesting feature of all TWs found so far is that they are unstable but 
only with an unstable manifold of very small dimension (invariably less than 10 
for those checked so far - Faisst & Eckhardt 2003, Kerswell & Tutty 2007). The 
typical situation is illustrated in figure 2 which displays the spectrum of the C3 
wave at the wavenumber a* which gives its lowest saddle node bifurcation. While 
Re is double- valued near the bifurcation, the phase speed C is not, monotonically 
decreasing in value as the bifurcation point is crossed from lower to upper branch. 
This then provides a convenient abscissa to show how the stability changes on both 
upper and lower branches as Re increases away from the saddle node bifurcation. 
Only disturbances which share the R2-symmetry of C3 are considered. C3 is par- 
ticularly interesting for two reasons: 1) the number of unstable directions decreases 
down to one for Re > 1826 along the lower branch and 2) there are several bifur- 
cations involving real eigenvalues. The implication of the first observation is that 
C3, which is on the laminar-turbulent boundary (in a pipe 2.5D long at Re = 2400, 
Duguet et al 2008a), will become a local attractor there beyond Re = 1826 for 
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Re 

(upper branch) 

1456 1195 



1146 



Re ► 

(lower branch) 

1261 



1575 



2015 



0.03 



0.02 



^ 0.01 



0.00 



-0.01 




1.36 



1.40 



1.44 



1.48 



1.52 



Figure 2. Stability of C3 within the R2-symmetric subspace against the phase speed C 
as Re increases away from the saddle node bifurcation at 1141 (q = a* and the corre- 
sponding Re is across the upper x-axis). Each line either indicates the locus of a real 
eigenvalue (solid) or a complex conjugate pair (dashed) as Re changes. The blue (dark) 
lines correspond to those which are symmetric under S while the green (light) lines are 
anti-symmetric under S. The vertical dotted line indicates the saddle node, of which the 
inset is a close up. (The dotted green indicates a stable complex conjugate pair which was 
difficult to resolve). 

R2-symnietric flow (confirmed by Duguet et al. 2008a). The second observation 
means that new solutions bifurcate from C3 which are also steady in an appropri- 
ately translating Galilean frame - that is to say they are traveUing waves (this is 
how the asymmetric travelling wave of Pringle & Kerswell 2007 appears). There 
are 4 possibihties: a transcritical bifurcation where no symmetry is broken, and 
3 types of symmetry-breaking pitchfork in which only S-, Z- or J72-symmetry is 
retained (recall 2 symmetries imply the third). For example, the bifurcation at 
{C,Re) = (1.47, 1449) is a mirror-symmetry-breaking pitchfork which can be fol- 
lowed to produce TWs which resemble those of S2. D2 is plausibly the product of 
a S-symmetry-breaking pitchfork from C3 (D2 is very similar to Z2 and hence C3 
- see figure 1). 

Hopf bifurcations are the generic scenario, however, leading to more compli- 
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Source 



Symmetries Reiou 



FE03 & WK04 



WK04 



PK07 



here 



S2 


S &R2 


1358 


1.55 


S3 


S &R3 


1251 


2.44 


S4 


S & R4 


1647 


3.23 


S5 


S & Rs 


2485 


4.11 


S6 


S &R6 


2869 


4.73 


? 


S & fii 


3046 


2.17 


Asymm (SI) 


S 


^820 


^ 1.8 


Mirror- sy mm (Ml) 


S & fii 


773 


1.44 


Helical 


n{ 


773 


1.44 


M2 (A3) 


S, R2 & ^2 


1125 


2.0 


M3 


S, Iv3 & ais 


1552 


2.2 


M4 


S, 1x4 & ^^4 


1824 


2.6 


M5 


S, Rs & fls 


2143 


3.1 


M6 


S, Re & fie 


2531 


3.5 


N2 (C3) 


S, R2 & fl2 


1141 


1.2 


N3 


S, R3 & ^3 


1037 


2.0 


N4 


S, R4 & O4 


1290 


2.5 


N5 


S, Rs & fls 


1622 


2.9 


Z2 


Z & R2 


^ 3250 


0.8 


D2 


Z & R2 


<2875 


? 



Table 1. The various symetries of all the TWs currently known. FE04 is Faisst & Eckhardt 
(2003), WK04 is Wedin & Kerswell (2004) and PK07 is Pringle & Kerswell (2007). The 
'?' mark in the entry for WK04 is to indicate that they found a mirror-symmetric TW 
but at high Re and its relationship to Ml is unclear. 

cated, relative periodic orbits such as that traced by Duguet et al. (20086). Within 
the Re range of figure 2, C3 experiences 3 Hopf bifurcations which all lead to relative 
periodic orbits with the Z-synimetry broken. 



4. New mirror-symmetric classes 

The fact that A3 and C3 appear at such low Re strongly suggests that there are 
analogous waves in different rotational symmetry classes also existent at pre- and 
transitional Re < 2400. This indeed turns out to be the case with all except one 
(M3) of the new families of mirror-symmetric modes being easily found using the 
homotopy approach once the appropriate ri-symmetry is incorporated. Selecting a 
roll forcing of R2m-symmetry defined, in the notation of Wedin & Kerswell (2004), 
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Figure 3. Slices of the lower branch M-class TW solutions at Re = 2400. M2 (A3,top 
left), M3 (top right), M4 (bottom left) and M5 (bottom right). Contour levels of the 
streamwise velocity perturbation are in increments of 0.02L*' running from — 0.19[/(dark 
red) to 0.19C/(white). 



by A2mi with i = 2,3 or 4 invariably led to a subharmonic Rm-symmetric streak 
instability which could be easily tracked back to zero forcing (this was exactly the 
strategy used by Wedin & Kerswell to find their Ri TW which has two roll pairs 
and appears beyond Re = 3046). The remaining wave M3 was found by homotopy 
in m starting at a M5 wave (M5 being used rather than M4 because of the similar 
parity of velocity components with respect to the radius). 

Inspection of the new lower branch TWs - see Figures 3 and 4| - indicates 
that the mirror-symmetric family reported in Pringle & Kerswell (2007) is the first 
family (Ml) and A3 is a member of the second family (M2) of a class of TWs (Mm 
with m=2,3,4,5,6, . . . ) with two layers of fast and slow streaks. Furthermore, C3 
appears a member of a R2-symmetric family (N2) of another class of TWs (Ntti 
with 771=2,3,4,5,...). The ordering of the families within the respective classes is 
clear from the phase speed and wavenumber data shown in figure 5. 



t Videos to be included with this paper are particularly illuminating 

http: //www. maths. bris.ac.uk/'^cpl571/TWs/table. htm 



they are currently on 
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Figure 4. Slices of the lower branch N-class TW solutions at Re = 2400. N2 (C3, top left), 
N3 (top right), N4 (bottom left) and N5 (bottom right). Contour levels of the stream- 
wise velocity perturbation are in increments of 0.02U running from — 0.19C/(dark red) to 
0.19(7(white). 

Extending the Mm and Nrn families to higher rotational symmetry m is straight- 
forward and could have been continued if the general trend had not already emerged. 
Finding lower families is, however, more difficult with Nl noticeably absent at the 
time of writing. This family could potentially be the first to appear as Re increases 
from zero but this seems less likely after plotting the optimal slices of each family's 
solution surface together on a phase speed verses Re plot ('optimal' being defined 
by that wavenumber a* which corresponds to the lowest saddle node bifurcation). 
Whereas the bifurcation points for Mm TWs monotonically increase in Re and 
decrease in C with m (see figure 6), this is not true for the N-class. Of N2 to 
N5, N3 actually has the lowest bifurcation point. This mimics the situation for the 
original S-class (non-mirror symmetric) TWs. The fact that the lowest Re for each 
Nrn family as well as Ml and M2 is consistently smaller than the lowest Re for 
the equivalent Sm family is suggestive that the latter all bifurcate off the former 
in mirror-symmetry-breaking bifurcations. This, after all, is the situation for the 
asymmetric waves reported in Pringle & Kerswell (2007) which form the missing 
(trivially) Ri-symmetric family from the original studies and which bifurcate off 
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Figure 5. A plot of the phase speed C (in units of L^) against the axial wavenumber a (in 
units of 2/D). The contours correspond to different Re with the loops (M-class/N-class 
TWs shown using blue solid/red dashed lines) constricting as they move towards the 
saddle nose indicating a* . The various Re for each TW are as follows: Ml 776 & 820, M2 
1125 & 1145, M3 1552 & 1650, M4 1846 & 2037, M5 2148 & 2400, M6 2540 & 2662; N2 
1150 & 1318, N3 1038, 1050 & 1120, N4 1292 & 1300, N5 1629 & 1652. 

the Ml family. For a given {a, Re), a TW needs at least 2 quantities such as the 
the phase speed C and a Reynolds number, Rep, measuring the wall shear stress to 
characterise it. Examining N3 and S3 TWs shows that they occur in similar parts 
of {Re, C, Rep) space at a = 2.44 (a* for S3) but there is no obvious connection 
between them at least for Re < 2500. The implication then is, if S3 does bifurcation 
off N3, the solutions exist to significantly subcritical Re which would explain why 
the importance of mirror-symmetric TWs hadn't been realised before. 

The new upper branch TW solutions - see figure 7 - all show an expected inten- 
sification of the slow and fast streaks, and the positioning of the fast streaks closer 
to the pipe wall than the corresponding lower branch solutions. This localisation 
is particularly noticeable for the N-class TWs indicating high wall shear stresses 
which is borne out by figure 8. The friction factors achieved by the upper branch 
and even by some of the lower branch N-class TWs as well as Ml and M2 are sig- 
nificantly higher than time-averaged experimental values (as shown by the upper 
dotted line) and those of the S-class TWs. 

Finally, it's worth remarking that the M-class TWs start to appear at fascinat- 
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Figure 6. A plot of C against Re for a* (see Table 1) which gives the lowest saddle node 
bifurcation for each TW (M-class/N-cIass waves shown using blue solid/red dashed lines). 
The non-mirror symmetric S-class of Faisst & Eckhardt(2003) and Wedin & Kerswell 
(2004) which generally appear at higher Re are shown using black dotted lines. Typical 
resolutions used (in the nomenclature of Wedin & Kerswell (2004)) are: Ml (14,25,5), 
M2 (12,30,6), MS (10,30,10), M4 (10,30,8), M5 & M6 (10,35,8); N2 (12,30,6), N3 & N4 
(10,30,8), N5 (8,35,8). 

ingly low Re given their intricate double radial layer structure: for example, M4 
appears at Re — 1824 (see figure 6 and Table 1). That they exist at all is not a 
surprise (presumably TWs with three radial layers are possible too) , but that they 
appear so early in Re surely is and contrasts starkly with Z2 which doesn't emerge 
until Re « 3250. 



5. Discussion 

In this paper we have described 2 new classes of mirror-symmetric TWs - the M- 
and N-classes which are also both shift-&-reflect symmetric - and a new mirror- 
symmetric family Z2 which is not. The M- and N-class waves appear earlier in 
Re than the original non-mirror-symmetric waves of Faisst & Eckhart (2003) and 
Wedin & Kerswell (2004) suggesting that the latter are borne through generally 
supercritical symmetry-breaking bifurcations from them. This was already found 
to be the case for the Ml waves (Pringle & Kerswell 2007) but now seems more 



Article submitted to Royal Society 



Highly- symmetric travelling waves in pipe flow 



13 




Figure 7. Slices of upper branch solutions at Re = 2400. M3 (top left), M4 (top right), 
N3 (bottom left) and N4 (bottom right). Contour levels as in figures 3 and 4. 

generally true now that the various TW families making up each class have been 
found. The stability analysis of N2 presented here provides a timely reminder of the 
bifurcation possibilities: S or Z symmetries can be broken at pitchfork bifurcations 
leading to either Z-class TWs or 5'-class TWs found originally (Faisst & Eckhardt 
2003, Wedin & Kerswell 2004). Hopf bifurcations, of course, give rise to periodic 
orbits in the Galilean frame of the TW or quasiperiodic orbits in the pipe (labo- 
ratory) frame. Tracing these requires a more sophisticated approach based upon 
time-stepping the equations and searching for an exact return of an appropriately 
chosen Poincare map (Viswanath 2007,2008, Duguet et al 2008b). 

The original motivation for searching for all the new TWs discussed here was 
the discovery of previously-unknown TWs in the laminar-turbulent boundary by 
Duguet et al (2008a) . It is likely that some of the lower branch solutions of these 
new waves similarly populate this boundary so that their stable manifolds also play 
a role in determining if a given initial condition will lead to relaminarisation or a 
turbulent evolution. There are, however, many interesting issues surrounding this 
assumption. For example, if a lower branch TW is on the boundary at one Re, will 
it still be at 10i?e and if not, how did it leave? In a long pipe where the TWs are 
continuously parametrised by their wavenumber over a finite range, at what critical 
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Figure 8. The friction factor, A, against Re for the newly found travelling waves, as 
well as the asymmetric wave (Ml). The lower dotted line corresponds to the laminar 
state, and the upper dotted line to the log-law parameterisation of experimental data, 
1/\/A = 2.01og(i?emVA) - 0.8 (see Schlichting 1968, equation (20.30)). The inset shows 
the same plot but for the travelling waves S2— S6 (only selected curves drawn for clarify) 
and Z2 which appears at much higher Re. The earlier onset in Re, and significantly higher 
friction factors of the new M and A^-class travelling waves is clearly apparent. 

wavenumber does the lower branch TW leave the boundary on its way to becoming 
an upper branch TW? The asymmetric wave (Pringle & Kerswell 2007) provides 
an obvious example being on the boundary for a = 0.625 (at Re — 2875, Schneider 
et al 2007b, Duguet et al 2008a) but presumably not for a = 1.44 where its wall 
shear stress is high (see figure 6 of Pringle & Kerswell 2007) . Hopefully these issues 
will be discussed elsewhere in this celebratory volume. 

The discovery that the upper branch solutions of the new N-class waves have 
such high wall shear stresses is, however, potentially more important. A recent 
attempt to extract coherent fast-streak structures from pipe turbulence (Willis & 
Kerswell 2008) has concluded that the TWs currently known were not energetic 
enough to be part of the turbulent attractor as previously supposed and speculated 
about others as- yet-undiscovered. The new N-class of mirror-symmetric TWs may 
well be these missing waves and the issue clearly needs to be revisited. 

In this report, we hope a step forward has been made in appreciating the hi- 
erarchy of TWs which exist in pipe flow. Generally the picture is that TWs with 
shift-&-reflect symmetry, mirror-symmetry and a low degree of rotational symmetry 
seem to appear first and spawn further, less symmetric, TWs through bifurcations 
as Re increases. However, given the notorious complexity of the Navier-Stokes equa- 
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tions, it would be foolhardy to be too rigid about this especially as the Nl family 
still remains 'at large'. What should be absolutely clear, though, is that the pipe 
flow problem continues to fascinate and intrigue us fully 125 years after Reynolds' 
original experiments. 

CP acknowledges studentship support from EPSRC and YD the award of a Marie-Curie 
IntraEuropean Fellowship (grant number MEIF-CT-2006-024627). 
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